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NATIONAL ADVISORY COMMUTES FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1280 


(21APEICAL DEIERMINAIION OF WALL TEMPERATURES 
FOR HEAT TRANSFERS THROUGH WAITS 
OF ARBITRARY SHAPE* 

By Otto Lutz 


A graphical method is given which permits determining of the 
temperature distribution during heat transfer in arbitrarily shaped 
walls. Three examples show the application of the method. 

The further development of heat engines depends to a great extent 
on the control of the thermal stresses in the walls. The thermal 
stresses stem from the nonuniform, temperature distribution In heat 
transfer through walls which are., for structural reasons, of various 
thicknesses and sometimes complicated shape. Thus, it is important 
to know the temperature distribution in these structural parts. Follow- 
ing, a method is given which permits solution of this problem. 


STATEMENT OF THE PROBLEM 


According to figure 1, a two-dimensional heat flow through the wall 
is assumed; thus, the wall should extend in the direction normal to the 
flow sufficiently uniformly so that no components of the flow occur in 
this direction. Furthermore, we consider the steady state process only. 
It should be added that an approximately steady state flow is present 
even in reciprocating engines since the temperature oscillations at the 
wall surface are only small and are rapidly damped in the wall.-*- 

*”Zeichnerlsche Ermittlung der Wandtemperaturen beim Warmedurchgang 
durch Wande von beliebiger Form." Zeitschrift des Vereines deutscher 
Ingenieure, Band 79* Nr. 34* August 1935* PP- 1041—1044. 

-^Compare G-. Eichelberg: Temperaturverlauf und Warmespannungen in 

Yerbrennungsmotoren. VDI-Forschungshef t 263* Berlin 1923. 
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The problem is characterized by the fact that in the region limited 
by the two boundaries and 02 ( wall -surfaces) and the two heat 

streamlines sp and S2 (the course of which is assumed to be known) 
(fig. 2), Laplaces* s equation for the potential distribution 


V 2 d 


_ agg , afo 
= bf 


0 



is valid, with i3 denoting the temperatures. The heat flow enters 
through the boundary Op, and leaves through 02- However, these two 
boundaries are not isotherms, but on the contrary, the boundary condi- 
tion of the third kind 


“l ( S 1 “ 13 1) 


= - X 


a n 


02 (g2 - 6 2 ) 




( 1 ) 


applies tp them, wherein a denotes the heat— transfer coefficients, 0, 
the outside temperatures, -3, the wall temperatures, and X , the con- 

c ^ 

ductivlty of the wall; — is the temperature gradient normal to the 

8n 

surface (directed inward). 


So far, no method giving a general solution of this problem is 
known to me. For the present case of heat transfer, Geiger 2 , Eichelberg, 
(footnote l) and Lachmann3 have striven for the solution by assuming the 
flow pattern (potential flow) and varying it until it corresponded to the 
boundary conditions, equation (l) . This method is possible for simple 
forms only, and, if carried out correctly it is rather slow. 


We take the inverse procedure: Assuming the flow direction through 

the boundaries, we first satisfy the boundary conditions (l), and check 
whether the flow pattern corresponds to a potential flow. The advantage 
of this method is that it leads to the solution comparatively quickly 
and directly shows the surface temperatures, which are of greatest 
interest. 

- •— ' : 

J. Geiger, Z. VDI Bd. 67, 1923, p. 905. 

3 k. Lachmann, Z. VDI Bd. 72, 1928, p. 1127. 
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. DEVSL0PME1W Off THE METHOD 


If we assume the flow direction through the surfaces and require 
that the same heat quantity flow between any two flow lines thus started, 
we shall, in general, obtain the flow pattern shown in figure 3 a if the 
flow direction was not, by chance, assumed correct initially. The flow 
direction is to be varied so that corresponding flow lines run into one 
another. Since, in most technically important cases, the thickness of 
the wall through which heat transfer takes place is small as compared 
to the dimensions along the surface, one will almost always get by with- 
out a determination of the potential field in the wall . It is sufficient 
if the condition shown in figure 3"b is reached wherewith the problem is 
solved. 

Assume that the flow entering q(kcal/n£h) forms the angle cp 
(measured positive in clockwise direction) with the perpendicular 
(fig. k). ffor the flow components in the direction of the surface o, 
q^, and q n In the direction n normal to the surface the equations 


( 2 ) 


90 

- - & 


are generally valid. They are connected with each other by 


- = q Q tan cp 


(3) 


If we derive from the boundary' condition (l) the temperature gra- 
dient along the surface, to which corresponds the flow qQ, the latter 
becomes 


_ i dq n 

110 a Bo 


or' with equation (3) 
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For the performance of the graphical. solution, we turn from the 
differential form to the difference form and obtain finally, at the 
surface through which the flow enters (subscript 1 ) 




tan cp-^ 


( 5 ) 


at the surface through which the flow leaves (subscript 2 ), according 
to figure h 


= © % % 


(5a) 


Therewith, a step— by— step determination of the flow distribution along 
the surface is possible- It will be useful to select the steps Ao as 
integral parts of the "fictitious transfer wall thicknesses" l/a, thus 
for instance, as 1 / 5 , l/lO, or if 20. (According to equation (5), the 
fractions must be chosen smaller, the larger tan 9 , thus, the more 
oblique the angle of incidence between surface and permeating flow. ) 


We now develop the two surfaces (fig. 5 ), assume the flows q^ 

and qg at the beginning of the developed regions at arbitrary linear 

lengths, and divide the developed surface into equal parts Ao. Start- 
ing from the center A of these parts, we draw the lines ABC given 


by the assumed flow angle 9 -, 


Aon 


(at A) and the slope tan ^ = — 7 I— 

M x/ a i 


(at B). BD is a parallel to the base line at the distance 


V 


The slope tan ^ is constant for the remainder of the procedure, since 
equal parts Ao^ are chosen. According to equation ( 5 ), CE then rep- 
resents the increase (here the decrease) Aq^ 


which, in addition, is 

Ao^ considered (FG). The 
describes the oourse of the flow permeating the surface. 


to be shifted to the end of the increment 
line IF . 


As is well known, for flow fields free from sources> the same 
partial flow must be transported in each "stream tube" or "flow tube." 
Accordingly, we divide the area EFEEK, representing the total, flow, 
into t equal area parts (Q^/t) and thus obtain the points of pene- 
tration of the (t— 1) flow lines at the surface 0 ^. The direction of 
the flow lines at these points is given by the assumed orientation of 
the field. 
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Repeating the construction at the exit surface ‘Oq (lines A* B*C’, 
etc. ), we obtain a diagram similar to figure 3a. The field orientation 
at the surface is to he changed so that corresponding flow lines meet. 

We shall later on find a few criteria according to which the corrections 
may he suitably performed. 


Once the flow pattern is plotted in this manner* there remains to 
he determined what actual flow passes through for a given difference of 
the outside temperatures . So far. We hnow only the flow distribution 


along the surfaces and the ratio — of the entering and leaving flows 

lug 


of a flow tube resulting from the requirement that the total entering 
flow (area DHIK) must equal the leaving one (area D’H’I’K'). The ordi- 
nates q^ would have to be changed in the ratio of these areas if they 

are to be comparable to the corresponding ordinates . 


The wall is cut open along the canter flow line s m of an arbi- 
trary flow tube (fig. 6). Since, in general, the flow q, according to 


q = 



(2a) 


corresponds to the temperature gradient and, according to figure k, 
q = q-n/toB cp, we plot a gradient corresponding to this flow q n /cos cp 
at the two surfaces. The appropriate procedure, according to figure 6, 
will be to represent the gradient as ratio of the flow q n _^/cos cp-j_ 

found, according to figure 5, to the total flow % represented by DHIK. 
Therewith, the above-mentioned transformation of the flows 


and 




in the ratio of the areas Q,^ and Qg becomes necessary. 


The scale is to be selected so that one obtains slopes usable in tbe 
diagram (30° to V?°) . Again both gradients are joined corresponding 
to the flow pattern, wherein intermediate slopes must be chosen 
inversely proportional to the flow widths in the flow pattern. 


Now, one has at the surface o-^ 

a l ( 0 1 “ *]) = qjai = 1 cos qp-j^ 


(6) 


or 


Q 1 ~ «1 
b 003 *1 



(6a) 
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Accordingly, the lines Indicating the gradient are extended outside of 
the vail up to the distances cos 9^ and ■ cos cpgjL then the differ- 
ence in height h corresponds to the prescribed temperature difference 
©1 — ©2* As a proof, this construction will -be carried out on two 
different flow tubes. 


Therewith, the temperature $ - x * at the point Qq is known. 

Since, according to the boundary condition, equation (l) or equation (6), 
a linear relationship exists between the flow q through the surface 


and the surface temperature, the flow distribution found., in figure 5 
also reproduces the temperature distribution at the surface. The 
temperature scale follows from the two known points X for =• 



thus = ©j_, and Y for the surface temperature Just determined 
8^* at 0 -^. The corresponding method is to be applied at the 


surface 02. 



EXAMPLES 

Wall with Protuberant and Indented Corner (Fig. 7 ) 


We presuppose that the flows develop in the two corners, indepen- 
dent of each other, so that the flow lines there may be assumed as lineB 
of symmetry at the corners. Only the region between these two flow 
lines 'will be considered. 

As sumpt ions . — Wall thickness, 25 mm; investigated length, 125 mm; 
heat conductivity, X = 50 kcal/mh°C; heat— transfer coefficients: 
a,j_ = 5OO, 0£ = 1000 kcal/n?h°C. 

Solution .— Surface o^: assumed flow qj_ ==-$6 mm 

Aoj_ 

Aoq = 10 mm, thus — ^ = tan 

Surface Oq : assumed flow qg 5, 25 mm; in order to 

obtain approximately the same values- Aq as for .surface a value 

of qg smaller than, q^ is chosen in Inverse ratio of til© heat- 
transfer coefficients. 


10 X iq -3 = 'q_ 
"30/500 10 
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Aog = 10 nm, thus 




10 X 10~~ 3 = 1 

50/1000 5 


The two flow areas are divided. Into t = 10 parts. -Then the 
surface temperatures •$ and ■Q 2*3 which establish the temperature 

scale 3 are determined: Since the surfaces are parallel the construc- 

tion given In figure 6 is simplified by determining those two points 
0 -l and 0^ located opposite each other where the perpendicular flow 

q n entering and leaving is the same. According to the flow pattern 
plotted in figure 7* the temperatures at these points are = lj2.8 0 , 

$2* = 28. 6°3 if the entire temperature difference is assumed to be 
9 ± - 02 “ 100 °- 


* Besides 3 these two temperatures would occur also in case of heat 

transmission through a plane wall of the thickness b = 25 mm. 

\ Result . — Thus , one obtains for the surface o^ the excess temper- 

ature at the protuberant comer as 63 0 — H2..Q 0 = 20.2°, or 
20 2 

Jqq — Q = 35-3 percent of the temperature jump, from the outside to 

the two-dimensional wall, and the insufficient temperature at the in- 
dented comer is 42.8° — 30-1° = 12-7° or 22.2 percent. 

At the surface 02 , the excess temperature at the indented corner 
amounts to ^1.2° — 28.6° = 12-6° or ^ = kk- percent of the tempera- 
ture Jump, from the two-dimensional wall to the outside, and the insuf- 
ficient temperature at the protuberant corner 28.6° — 1^. 3° = 14-3° 
or 50 percent. The differences are relatively larger on the side with 
the larger heat— transfer coefficient. Moreover, it can be inferred 
that the flow lines turn their concave side more strongly toward this 
surface. 


Protuberant Comer Rounded Inside (Pig. 8) 

The same dimensions and material constants are chosen as in the 
first example. At a sufficient distance from the corner, the heat flow 
will pass through the wall at a right angle; there the entering flow 
then equals the leaving flow qg. If the graphical determination 
is started at this point with q^_ = q2> the two flow areas and Q2 

resulting, according to figure 5, must be equal. 
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Surface o-^: assumed flow =- 50 ram. r 


Aoj_ = 10 mm. 




10 X 10~3 = 1 
50/500 10 


Surface Og: assumed flow = 50 mm 


AOg = 10 mm.j 


thus 



10 X 10^3 1 

50/1000 5 


Temperature scales : $ and i3 g* are determined at the starting 

point of the investigation. 

R^ult.- Surface o 1 t Excess temperature at the salient corner 

^*9 _ ~ 42.8° = 24.1° or 42.1 percent of the temperature jump from 
outside to the two-dimensional wall. 

Surface o 2 : Excess temperature 36.8° - 28.6° = 8.2° 

or 28.6 percent. . * 


Due to the rounding off, the temperature differences become larger 
outside, and smaller inside. 

Wavy Wall (Jig. 9) 

The waviness is assumed as sinusoidal; 1 = 50 kcal/mh°C: a n = -575 
Og = 1250 kcal/m 2 h°C. ’ ■ - 1 - 


Surface o-j_: assumed flow = 20 


mm 


*°i ’ ^ (• t i) - " °- 03 


Surface og: assumed flow qg = 20 


mm . 


Acg = 4 mm, thus — ft (= tan fg] = — X , 10 — = 1 — 
A/ctg \ / 50/1250 10 
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Temperature scales: The flow pattern for determination of S • 

and -Sg* according to figure 6 is plotted for the "boundary stream, 
lines . Sj_ and Sq. 

Result .— The temperature differences are only minor at the surface 
with the smaller heat— transfer coefficient, "but' considerable at the wavy 
surface where., moreover, the larger heat— transfer coefficient prevails. 

From the graphical construction and the examples , the following 
criteria for correction of the assumed flow direction may be derived: 

1. The concave side of .the flow lines is turned more strongly 
toward the surface with the larger heat— transfer coefficient, or 
the surface with- the larger heat— transfer . coefficient is crossed 
by the flow lines at a steeper slope providing other conditions are 
equal. 

2. The width of the flow tubes along the entrance surface de- 
crease in the direction toward which the flow lines slope. (Coup are 
figs. 7 and 8 . ) 

3 . The width of the flow tubes at the exit surface increases 
toward the side in which the flow direction is sloping. 

1. The differences in the widths along a surface are greater, 
the more oblique the angle at which the flow lines cross the surface. 

5 . Other conditions being equal, the differences in the widths 
along fhe surface are greater, the larger the heat transfer 
coeff icient.^ 


Translated by Mary L. Mahler 
Rational Advisory Committee 
for Aeronautics 


^Forschg. Ing.— Wes. Bd. . 6 , Rr. 5, 1935, p. 2kO } reports on an extension 
of the method to include axi ally symmetrical forms, locally different heat- 
transfer coefficients, and locally different outside temperatures. 
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Figure 1.- Two-dimensional heat flow through a wall. 



Figure 2.- Section limited by the wall boundaries 0 } and 03 and the heat 
streamlines and Sg. . The heat flow Q, enters through the 
boundary o^ and leaves through the boundary 0 £. 
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Figure 3.- Heat -flow patterns in case of incorrect (a) and correct (b) 
assumption of the flow direction through the surfaces o^ and o 2 - 



Figure 4.- Decomposition of the entering and leaving flow into its 
components in Hie direction of the surface and normal to it. 
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Figure 5.- Graphical determination of the course of the entering and 
leaving heat flow over the developed surfaces . 
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Figure 7.- Heat transmission through a wall with protuberant and indented 
corners. Temperature difference between the two sides e-^ - 02 = 100°; 

surface temperature on the entrance side, $ 2 on the exit side; 
heat conductivity A = 50 kcal/mh°C; heat-transfer coefficients = 500, 
02 = 1000 kcal/m2h°C. 



Figure 8.- Heat transmission through a protuberant corner with rounding oii 
the inside. Outer temperatures and material constants as in figure 7. 
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Figure 9.- Heat transmission through a wavy waU. Difference of the outside 
temperatures e-j_ - 9g = 100°; heat conductivity A = 50 kcaI/mh°C; heat- 

transfer coefficients 0 ^ = 375, = 1250 kcal/m2h°C. 


NACA- Langley - 4-10-80 -050 







